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Evidence for a pseudogap in underdoped Bi2Sr2CaCu208+5 and YBa2Cu3O6.50 from 

in-plane optical conductivity measurements 
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(Dated: April 14, 2008) 

The real part of the in-plane optical self-energy data in underdoped Bi2Sr2CaCu208+4 (Bi-2212) 
and ortho II YBa2Cu306.5 contains new and important information on the pseudogap. Using 
a theoretical model approach we find that the density of state lost below the pseudogap Apg is 
accompanied with a pileup just above this energy. The pileup along with a sharp mode in the 
bosonic spectral function leads to an unusually rapid increase in the optical scattering rate and a 
characteristically sloped peak in the real part of the optical self-energy. These features are not found 
in optimally doped and overdoped samples and represent a clearest signature so far of the opening 
of a pseudogap in the in-plane optical conductivity. 

PACS numbers: 74.25.Gz, 74.62. Dh, 74.72. Hs 
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A striking feature of the underdoped cuprates, not 
present in conventional metals, is the appearance of a 
pseudogap [ij below a characteristic temperature T*. In 
general, the opening of a pseudogap Apg leads to a reduc- 
tion in the electronic density of states (DOS) around the 
Fermi energy. Since, as yet, there is no full understanding 
of the exact origin of the pseudogap in the cuprates, an 
important issue is how it should be modelled, its width 
in frequency, its behavior as a function of doping and 
temperature and what happens to the spectral weight in 
the gap. While the first evidence of a pseudogap in the 
cuprates came from NMR experiments [2, y] that showed 
a gap in the spin excitation spectrum, a gap in the charge 
excitations was found in the c-axis conductivity of un- 
derdoped YBCO by Homes et aL[J]. The conductivity 
in this direction is incoherent, and like c-axis tunneling, 
is expected to be proportional to the density of states. 
Homes et al. concluded that the spectral weight lost in 
the pseudogap region was recovered at much higher fre- 
quencies. A similar conclusion was reached by Loram et 
at from specific heat measurements [5|. In contrast, Yu 
et al. from a recent study of c-axis conductivity of under- 
doped RBa2Cu306.5 (R=Y, Nd, and La) extending the 
measurements to higher frequency found that spectral 
weight was recovered [gj]. Also tunneling data suggests 
that the "coherence peak" that signals a tunneling con- 
ductance that is enhanced over the background persists 
well into the normal state in Bi-2212 [7| suggesting a re- 
covery of gapped states. 

The pseudogap is more difficult to identify in the ab- 
plane transport. While it was clear very early on that the 
ab-plane absorption showed an unmistakable decrease at 
energies below Apg [8, 4 SO, (ill, llJI ; ^ clear interpretation 
in terms of a pseudogap was problematic since it is diffi- 
cult to separate the effects of a pseudogap (a pseudogap 
in the density of states) and a gap in the spectrum of 
inelastic excitations (a spin gap) since both reduce ab- 
sorption at low frequencies. Further, studies of the real 



part of the conductivity did not show the expected re- 
duction in optical spectral weight below the pseudogap 
energy [l3[. This remains unexplained. In this letter we 
find that by focussing not on features in the optical con- 
ductivity, but on the real and imaginary parts of the 
optical self-energy we are able to separate the contribu- 
tions of the pseudogap and inelastic scattering. Further- 
more, simulations with simple models show that spectral 
features in the density of states are less spread out in 
the self-energy than in the conductivity making small 
changes easier to estimate accurately. We can trace to 
missing states below the gap to the region just above it. 

It has become common in the optical literature to 
analyze reflectivity measurements in terms of the real 
and imaginary parts of the optical conductivity ct^lu) — 
ai{uj)+ia2{(^) and present the results for the optical scat- 
tering rate 1/t°p{uj) as a function of oj. The extended 
Drudc formula, which is valid for correlated electrons, 
can be written as ATTa{uj) = iu!p/[u! — 2T,°p{uj)], where ujp 
is the plasma frequency and S°P(a;) = 'S°^{uj)+iT,2^(uj) is 
by definition the optical self-energy. Its imaginary part 
gives the optical scattering rate and its real part is re- 
lated to the optical mass re-normalization. The real and 
imaginary parts of T,°p(lu) are related by Kramers-Kronig 
(K-K) transformations. They play, for the optical con- 
ductivity, a role similar to that played by the quasipar- 
ticle self-energy 'S'^p{uj) in angle-resolved photoemission 
spectroscopy but these are different quantities although 
they are closely related [iJ, |l5| . 

The real part of S°P(a;) obtained from the in-plane 
conductivity in the highly underdoped regime|12Lll6l.ll7| 
shows a particularly clear and striking signature of the 
pseudogap modelled just as was done from the specific 
heat [5] by a simple reduction in electronic density of 
states in the vicinity of the Fermi energy. Contrary to the 
specific heat case, however, we find evidence that these 
missing states in DOS pile up in the region just above 
the pseudogap energy. While the exact shape taken by 



the quasiparticle DOS around the Fermi energy will in 
principle be reflected in the shape of the real part of the 
optical self-energy, these details are of secondary impor- 
tance. 

In Fig. [T] we conipare results for — 2E°^(w) vs. uj ob- 
tained in ref. [16|, li2| for an underdoped Bi-2212 sam- 
ple (top frame) with an overdoped sample with approx- 
imately the same critical temperature (bottom frame). 
The evolution of the curves seen as the temperature is 
lowered is quite different in the two cases. The over- 
doped samples show only small changes while the un- 
derdoped one shows the growth of a sharp triangle-like 
peak around 750 cm~^. We will argue that this sharp 
peak is a result of a pseudogap opening with a recov- 
ery region just above it while at the same time a sharp 
resonance exists in the electron-boson spectral density 
I'^xi'^) which we believe to be due to spin fluctuations 
and not phonons. The inset in the top frame of Fig. 1 
shows results from the highly underdoped ortho II com- 
pound YBagCu-iOfi f^|12l| in which every second chain is 
full and the others are empty, leading to a perfect stoi- 
chiometric sample. The same characteristic temperature 
evolution is seen in YBCO as in Bi-2212. This suggests 
that this sharp peak is a universal feature of the under- 
doped cuprates. 

To establish our central point it will be sufficient to use 
approximate, but surprisingly accurate, formulas that re- 
late the real and imaginary parts of the scattering rates 
to the electron-boson spectral density Pxi'^) for zero 
temperature from the work of Mitrovic et al. 18| and 
others [ilHIIllliil. At T = we have^, 
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where N{lo) is the renormalized electronic density of 
states symmetrized about the Fermi energy. It contains 
the pseudogap and, for finite bands, would decay rapidly 
at the renormalized band edge. There is a cutoff fre- 
quency, f2c on the spectral density. Here we used Q.c — 
4000 cm-i . 

We begin with the following observation. For a con- 
stant density of states and coupling to a single mode 
of energy we i.e. Pxil^) = -^o^if^ ~ ^e) with ^o set- 
ting the scale, the optical scattering rate {1/t°p{ijj)) is 
zero for w < uje and equal to 2'kAq[{ijj — uje)/'^] for 
Lo > u>E which starts from zero at w = lue and satu- 
rates to 2ttAo only for uj 3> lue- If, in addition, there is 
a full gap Apg (the pseudogap) in the symmetrized den- 
sity of electronic states N{uj), 1/t°'''{uj) starts instead 
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FIG. 1: (color online) Minus twice the real part of the optical 
self-energy — 2E°''(a;) vs. uj for underdoped (top frame) and 
overdoped Bi-2212 (bottom frame) [Iq. [l7| . The inset shows 
the corresponding results for underdoped ortho II YBCO 13] • 



Apg and rise towards saturation 



UJE - ^pg)/'^ 



from zero at w = ue 

according to the modulating factor (uj 

which effectively rises more slowly with to than does 

[ij — uje)/lo. This is seen in the middle frame of Fig. 

2 where we show — 2E2^(a;) for three cases. A, B, C all 

based on the electron-boson spectral density shown in the 

inset, Px{^) = Ap/[\/2^(d/2.35)]e-('^-"^)'/[2(rf/2.35)=] ^ 

As wl\Jlf + ijj\ where A^ = 200, cjp = 250 cm"\ d = 
10 cm^^, As = 100 cm^^, and Ws/ — 500 cm^^. It has a 
peak at w = 250 cm~^ and a broad background extend- 
ing to 4000 cm~\ modelled on results for underdoped 
ortho II YBC0[13|- The curve A includes no pseudogap 
while curve B has Apg = 550 cni~^ and starts from zero 
at 800 cm~^ in contrast to curve A which starts at 250 
cm~^. Comparing with curve A shows that a pseudogap 
leads to an optical scattering rate which rises less steeply 
from zero than it would with Apg = 0. A slower rise 
would also arise if instead of using a delta function at 
uje for t^x^^) part of the weight (not shown here) was 
distributed to a region of higher frequencies. The two 
curves for —^Yf^iy) as well as C continue to rise at the 
largest photon energies shown because we have included 
a background in t^xil^) which extends to 4000 cm~^. 
The curve C includes a full pseudogap just as in curve B 
but now it is also assumed that the lost states in 7V(w) 
reappear just above Apg in the interval Apg to 2 Apg. 



A: No Gap 

B: Full Gap No Recovery 

C: Full Gap Full Recovery 




FIG. 2: (color online) Model results for -2E°''(a;) vs. uj (top 
frame) and —2'E2^{u!) vs. u; (middle frame). The electrons are 
coupled to bosons with electron-boson spectral weight shown 
in the inset of the middle frame. A, B and C curves are with 
no pseudogap, pseudogap without recovery region, and pseu- 
dogap with recovery, respectively. The corresponding density 
of states, N{lu), is shown in the inset (bottom frame) with the 
same line types. In the bottom frame we show the derivative, 
d[-2E°/{uj)]/du vs. CO. 



These extra electronic states lead to increased scattering 
and the curve C rises from zero much more rapidly than 
curve B. Also at the energy marking the end of the en- 
hanced density of states (recovery region) there is a kink 
in l/r°P(aj) beyond which the rate of increase in l/r°P((jj) 
becomes less rapid. 

In the top frame of Fig. [2] we show the real part of 
I]°p{llj) which corresponds to the K-K transformation of 
the optical scattering rates of the middle frame. For the 
curve A there is no singularity at uje- Rather there is a 
very broad maximum at w '--^ V^i^e followed by a slow 
drop to about 1/3 its maximum value at 4000 cm"^. If, 
on the other hand, we include a pseudogap without a 
recovery region above it, we get curve B, for which the 
peak falls at ~ v2(w_e + ^pg)- Including a recovery re- 
gion in the density of states changes the shape of E^^(ijj) 
significantly as we can see in curve C. Here we include 
a full gap below Apg and displaced states just above as 
shown in the solid (red) curve of the inset on the bot- 
tom panel of Fig. [2] where our model N{oj) is shown. In 
this cases — 2I]°^(w) develops a much more pronounced 
peak between uje + Ap„ and loe + 2Ap„, points marked 
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FIG. 3: (color online) Comparison of —2'E°^{uj) vs. u for 
underdoped Bi-2212 and overdoped Bi-2212 samples, bot- 
tom frame. The top frame shows results of model calcula- 
tions. Solid (dark-blue) and dashed (blue) curves are with 
pseudogap and recovery and without gap (see corresponding 
quadratic behavior N{uj) curves in the inset) for the I^xi^) 
shown in the inset to Fig. (2] Dashed (purple) and solid (green) 
curves are obtained when the electron-boson spectral density, 
just has a continuous background. 



by maximum slopes. This structure is seen both in un- 
derdoped Bi-2212 and YBCO samples at low tempera- 
ture normal state of Fig. [T] and is strong evidence for a 
pileup of the spectral weight lost below Apg, in the re- 
gion Apg < uj < 2 Apg. The above characteristics are a 
clear signature of pseudogap formation in the in-plane 
conductivity. 

The feature of the model curves displayed in the top 
frame of Fig. [5] which we have just described and which 
are most important for the interpretation of experiments, 
shows up even more clearly in the derivative of the real 
part of the optical self-energy d[~2'Si^ (uj)] / duj . Results 
are shown in the bottom panel of Fig. [2l The curve A 
shows maximum slope at the energy uje of the peak in 
/^x(ti;) (see inset in the middle frame of Fig. [5]) and it is 
zero at ~ \/2u!e. Curve B has a maximum at uje + Apg 
and a zero at ^ v2(w_E + Apg) while curve C peaks at 
UJE + Apg and shows a minimum at uje + 2Apg indicating 
the end of the recovery region in our model of the effective 
electronic density of states, N{uj). 

In the lower panel of Fig. [3] we repeat our experimen- 
tal results for the underdoped and overdoped samples of 
Fig. [T] but for clarity, show only the two lowest temper- 
atures. In the top frame we show four curves. One set 
is based on the electron-boson spectral density shown in 



the inset of Fig. [2l The other is obtained when only the 
background is 

kept in /^x(a;). In Bi-2212 series we have found that 
in the overdoped regime there is no prominent peak in 
the electron-boson spectral density, only a background 



M llM iM |2J]. For this figure we have used 
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quadratic pseudogap shown as the (dark-blue) solid curve 
in the inset of Fig. [3l Starting with the two curves for 
the background only, Pxi^) = 100a;/[500^ + w^], we see 
that including a pseudogap does not change the curve for 
S^''(w) much in agreement with the data and there is, in 
this sense, no clear signature of Apg in such spectra. On 
the other hand contrasting with the other two curves, the 
pseudogap now clearly changes the shape of Ej^(aj) and 
brings it much closer to the experimental results for the 
underdoped sample (lower panel) . The peak in the (dark- 
blue) solid curve can be made sharper and resemble even 
more a logarithmic singularity by decreasing the width 
of the recovery region. This arises because the sharp rise 
seen in 1/t°p{uj) (curve B) of Fig. [2] middle frame, can 
be made even steeper and come closer to a pure vertical 
rise. We conclude that the striking difference observed in 
T,°^(uj) between underdoped and overdoped samples can 
be understood to result from the absence of a pseudogap 
and a sharp peak in I^xi^) for the overdoped case, and a 
combination of a peak and pseudogap in the underdoped 
case. Both are needed to get agreement with the data, 
including a recovery of spectral weight in the electronic 
density of states in the region just above the pseudogap 
with the end of the recovery region marked by a kink 
in E°^((jj) around 940 cm~^ for Bi-2212 (purple arrow in 
lower panel of Fig. [3]) and similarly 970 cm~^ for YBCO. 
If we assume the position of the peak in I'^xi'^) to be set 
by the neutron resonance energy given by fires = ^AksTc 
[25| we conclude that the pseudogap has a value of 43 
meV in Bi-2212 and nearly the same 45 meV in YBCO 
and that the recovery region is approximately Apg wide. 
We note that our suggestion that the recovery of spectral 
weight occurs immediately above the gap region in the 
normal state just above Tc, is in accord with tunneling 
spectra \J\. 

In summary, the real part of the optical self-energy, 
— 2E°^(w), defined through a generalized Drude formula 
at low temperature in the normal state, for the un- 
derdoped case, shows characteristic triangular behavior 
around 750 cm~^ not seen in overdoped samples. This 
behavior is traced to the existence of a pseudogap Apg 
in the electronic density of states with states lost below 
Apg recovered in the region immediately above it as well 
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as the existence of a peak in the electron-boson spectral 
density responsible for the inelastic scattering. By con- 
trast, in the overdoped samples no pseudogap opens and 
the fluctuation spectrum /^x(tj) has no resonance peak 
but rather consists only of a reasonably flat background. 
Criteria to recover a value of the pseudogap and the end 



of the recovery region are given and applied to under- 
doped Bi-2212 and ortho II YBCO. The signature of the 
pseudogap is much clearer in the optical self-energy than 
it is in the in-plane conductivity, the quantity used in 
previous studies. 
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